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In this proceedings paper, I overview the current status of the X-ray luminosity 
function of AGNs in the soft (0.5-2 keV) band, extended using XMM-Newton and 
Chandra survey data. We found that the number density of low luminosity AGNs 
peaks later in the history of the universe (z ~ 1) than the that of high luminosity 
AGNs {z ^ 1.7 — 3). I also describe the basic results of a spectroscopic foUowup 
fH ^ project of a complete HEAO-1 hard X-ray limited sample of AGNs using ASCA and 

^ ^1 XMM-Newton and present separate intrinsic hard X-ray luminosity functions for 

I ' unabsorbed and absorbed AGNs. We found that the absorbed AGN XLF drops 

O ' more rapidly at high luminosities, indicating a deficiency of absorbed luminous 

^ ! AGNs. 

1. Introduction 

' X-ray surveys are practically the most efficient means of finding active 

■ galactic nuclei (AGNs) over a wide range of luminosity and redshift. In 

order to construct an X-ray luminosity function (XLF) of AGNs, enormous 
efforts have been made to follow up X-ray sources with optical telescopes 
to establish their nature as AGNs and to measure their redshifts. Now that 
we have fairly complete samples of X-ray selected AGNs over 6 orders of 
magnitude in flux, from surveys ranging from all the high galactic latitude 
sky to the deepest pencil-beam fields. These enable us to construct and 
probe luminosity functions over cosmological timescales. 

In this proceedings article, I overview the current progress of a few 
projects related to AGN XLF. Firstly, I overview the results of the soft 
X-ray (0.5-2 keV) luminosity function (SXLF), which is the continuation 
of our previous work with iZ 05"^ T" samples^ ^, extended with deep XMM- 
Newton and Chandra surveys. 

While the soft X-ray surveys select against obscured AGNs, in the cur- 
rent situation, the number of available objects and area-flux coverage from 
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extensive surveys make them useful for probing detailed behaviors of XLFs 
of the unabsorbed portion of the AGN activity. A complementary XLF in 
the hard band (2-10 keV) (HXLF) enables us to also look into obscured 
AGNs, and thus it provides most direct measure of the accretion onto su- 
permassive blackholes (SMBHs). Ueda et al. in this volume covers our 
extensive recent work on HXLF (see also Ueda at al. 2003^° for a full 
description). 

In this article, we also present the results of our XMM-Newton and 
ASCA spectroscopic foUowup of a complete hard X-ray flux- limited sample 
of bright AGNs selected from HEAO-1 catalogs^^, the basic results of which 
have been integrated in the Ueda et al's HXLF. 

2. AGN Soft X-ray Luminosity Function and Evolution 

2.1. The Combined Sample 

In addition to the i?OS'^T samples used in our previous work^^'^^, we have 
added AGNs from the ROSATmTth Ecliptic Pole Survey (NEPS)^, from 
an XMM-Newton observation on the Lockman Hole ® and the Chandra 
Deep Surveys South^^/North^. For the Lockman Hole region, the in- 
ner part based on the ROSAT HRI has been replaced by a new XMM- 
Newton sample^. Four medium-deep 05^ T surveys used in our previous 
work^^'^^, i.e., the UK Deep Survey^, Marano Field^^, North Ecliptic Pole 
(deep PSPC-pointing)^ and the outer part of the Lockman Hole (PSPC)^^ 
are now collectively called the i?(95'^T Medium-Sensitivity Survey (RMS). 
The combined area versus limiting flux relation and the redshift-luminosity 
diagram are shown in Fig. 1. The samples are summarized in Table 1. 
We have tried to limit our analysis to "type 1" AGNs (including narrow- 
line Seyfert 1 galaxies). In the i? T samples, we selected type 1 AGNs 
mainly from the optical classifications. For the XMM-Newton Lockman 
Hole and Chandra DEEP field samples, the optical classification is supple- 
mented by hardness ratios of the X-ray sources. AGNs with z < 0.015 
have been excluded from the analysis to avoid the possible effects of the 
local large scale structure. Details will be explained in Hasinger et al. (in 
preparation). 

2.2. SXLF and Evolution with Redshift 

The soft X-ray Luminosity functions (SXLF) in different redshift bins have 
been calculated using the N°^^/N'^'^^ estimator, where each redshift bin 
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Figure 1. Left: The survey area of the combined soft X-ray sample as a function of 
fiux. Right: The soft X-ray sample in the z-Log Lx plane. 



Table 1. The Soft X-ray Sample 



Survey 


Area 


'S'xl4,lim 


Nagn 


Surey 


Area 


'5'xl4,lim 


A''agn 




[deg2] 


[cgs] 






[deg2] 


[cgs] 




RBS^« 


2 ■ 10'' 


250 


203 




1.-0.5 


.74-. 32 


83 


SA-Ni 


684.-35. 


47.-13. 


134 


LH-XMM** 


.33-. 10 


.13-. 08 


42 


NEPS*^ 


80.-0.35 


80.-1.5 


162 


CDF-S^^ 


.09-. 02 


.06-. 02 


115 


RIXOSio 


19.-16. 


8.5-3.0 


196 


CDF-N2 


.09-. 02 


.06-. 005 


97 



has been maximum-likelihood fitted with a smoothed two power-law form 
and the model value at the center of each bin is multiplied by the ratio of 
the actual number of AGNs in the bin to the model-predicted number^^. 
Nominally, corrections for incompleteness due to unidentified X-ray sources 
have been made by using an effective survey area, derived by multiplying 
the geometrical survey area by the completeness of the survey (i.e. iden- 
tified fraction of the detected X-ray sources). This method is valid when 
sources remain unidentified because of random reasons that are not corre- 
lated with the intrinsic properties of the source (e.g. optical magnitude). 
This is not necessarily true, especially in the deepest surveys. Thus we also 
calculated the XLF or number density upper bounds, where all the uniden- 
tified XMM-Newton and Chandra sources are assigned (in duplicate) the 
central redshift of each bin. Figure 2 (left) shows the SXLF in different 
redshift bins (plotted only the nominal incompleteness correction case). 
Figure 2 (right) shows the AGN number densities as a function of red- 
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Figure 2. Le/t; The SXLF for different redsfiift bins, calculated using the Afobs/Armdl 
method. Errors bars corresponds to Poisson Icr errors. Right: The number densities of 
soft X-ray selected AGNs with luminosities below and above Log Lx =44.5 ergs~^. The 
same curve for optically-selected QSOs with Mb < —26.0 are shown. The (S7m,f2A) = 
(1,0) cosmology is used in this plot for historical comparisons. 



shift separately for low and high luminosity AGNs, The incompleteness 
upper bounds (see above) are shown in dotted lines. Even in the most 
extreme cases of the incompleteness correction, we see that the number 
density of the low-luminosity AGNs peaks much later in the history of the 
universe (z ~ 1) than the high-luminosity case. This is in the opposite 
sense to the prediction from a analytical model based on the hierarchal 
merging and self-regulated accretion by Wyithe & Loeb According to 
their prediction, the number density of more luminous AGNs peak at later 
in the history of the universe. On the other hand, a prediction from a nu- 
merical simulation by Di Matteo et al,**, where gas density, star formation, 
and AGN formations are assumed to be related in a certain simple way, is 
consistent with the observed trend. 

For comparison, we overplot the redshift evolution of optically-selected 
luminous (Mb < —26) QSO number density^'''^. In the high X-ray lumi- 
nosity bin, we are not still certain (within the uncertainties in the incom- 
pleteness correction) whether we have detected the decline (with z) in the 
number density at z > 2.7, where the densities of luminous optical ^^'^ and 
radio^^ QSOs clearly show a drop. Note, however, that when we take a 
different method on incompleteness correction involving optical magnitude 
limits with our soft X-ray sample, a density decline at z > 2.7 is preferred 
even for the high-luminosity sample (Hasinger et al. in preparation). 
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3. The Brightest Hard X-ray Sample 
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Figure 3. Left: The local de-ahsorhedWKL¥ (2-10 keV) of unabsorbed (Log TVh cm-^ < 
21.5; filled circle) and absorbed (> 21.5; filled triangle) AGNs. The open squares show 
the sum of both. Right: Intrinsic absorption of the sample AGNs are plotted against the 
intrinsic (de-absorbed) luminosity. There would be ~ 9 objects in the region enclosed by 
a thick dashed line if absorbed and unabsorbed AGNs intrinsic HXLFs were the same 
besides normalizations. 



In order to construct the luminosity function of AGNs in the intrin- 
sic X-ray luminosity, X-ray surveys in the hard band [E > 2 keV) and 
X-ray spectroscopic (or at least hardness) information are crucial. In or- 
der to complement deeper ASCA^ XMM-Newton and Chandra surveys with 
smaller survey areas, we have defined a bright hard X-ray selected sample 
from the HEAO-1 all-sky surveys. In addition to the famous Piccinotti et 
al.^^ sample from the HEAO-1 A2 experiment, we have defined a somewhat 
deeper hard X-ray flux-limited sample of AGNs from the MC-LASS cata- 
log of X-ray sources from the HEAO-1 A1/A3 experiments in a limited 
region. The AGN sample from the latter catalog was investigated in detail 
by Grossan^. As a total, 49 AGNs are defined and we have made spectral 
analysis of all of them (except one) using ASCA and XMM-Newton obser- 
vation from archive as well as our own proposals. We have determined the 
intrinsic absorption A^h and underlying power-law index T. This enabled 
us to construct separate local HXLFs for absorbed and unabsorbed AGNs 
as functions of de-absorbed (intrinsic) luminosity (Fig. 3, left). We see that 



''http: / /heasarc. gsfc.nasa.gov/docs/heaol/archivc/hcaol_catalog. html 
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the absorbed AGN HXLF drops more rapidly than the unabsorbcd one at 
high luminosities. This can also be demonstrated in the - Nu plot (Fig. 
3, right). Suppose absorbed and unabsorbed AGNs had the same HXLF 
shape, there would be 9 AGNs in the region enclosed by a thick dashed 
line in this figure. See Shinozaki et al.^^. Full results will be reported by 
Shinozaki et al. (in preparation). 
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